Flow analysis with cumulants: direct calculations 
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Anisotropic flow measurements in heavy-ion collisions provide important information on the prop- 
erties of hot and dense matter. These measurements are based on analysis of azimuthal correlations 
and might be biased by contributions from correlations that are not related to the initial geometry, 
so called non-flow. To improve anisotropic flow measurements advanced methods based on multi- 
particle correlations (cumulants) have been developed to suppress non-flow contribution. These 
multi-particle correlations can be calculated by looping over all possible multiplets, however this 
quickly becomes prohibitively CPU intensive. Therefore, the most used technique for cumulant cal- 
culations is based on generating functions. This method involves approximations, and has its own 
biases, which complicates the interpretation of the results. In this paper we present a new exact 
method for direct calculations of multi-particle cumulants using moments of the flow vectors. 

PACS numbers: 25.75.Ld, 25.75.Gz, 05.70.Fh 



I. INTRODUCTION 

Anisotropic flow is a response of the system created 
in a heavy-ion collision to the anisotropics in the initial 
geometry. Thus, anisotropic flow is very sensitive to the 
properties of the system at an early time of its evolution. 
The sizable azimuthal momentum-space anisotropy ob- 
served at RHIC energies (for a review, see [HE]) is the 
main evidence for the nearly perfect liquid behavior [3J [3] 
of the created matter. Quantitatively, anisotropic flow is 
characterized by coefficients in the Fourier expansion of 
the azimuthal dependence of the invariant yield of parti- 
cles relative to the reaction plane [5[ [5] : 
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Here E is the energy of particle, pt is the transverse mo- 
mentum, (j) is its azimuthal angle, y is the rapidity, and 
^r the reaction plane angle (see Fig [l]) . The first co- 
efficient, Vi, is usually called directed flow, and the sec- 
ond coefficient, vi, is called elliptic flow. In general the 
v n = (cos[n((f> — ^rp)]) coefficients are p t and y depen- 
dent - in this context we refer to them as differential flow. 
The integrated flow is defined as a weighted average with 
the invariant distribution used as a weight: 
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Since the reaction plane ^ r is not known experimen- 
tally, the anisotropic flow is estimated using azimuthal 
correlations between the observed particles. For exam- 
ple, using 2-particle azimuthal correlations: 

<coB(n(& - fa))) = (e in ^~^) = (vl) + S n , (3) 



FIG. 1. Schematic view of a non-central nucleus-nucleus 
collision in the transverse plane. 



where the first term, (u 2 ), is the part due to anisotropic 
flow, and 5 n represents the so called non-flow contribu- 
tion, that comes from correlations not related to the 
initial system geometry. If non-flow is small, Eq. ^ 
can be used to measure v n , but in general the non-flow 
contribution is not negligible. To suppress non-flow one 
can exploit the collective nature of anisotropic flow using 
multi-particle correlations. The method based on multi- 
particle cumulants (genuine multi-particle correlations) 
to measure anisotropic flow was proposed in [7HTU] . This 
method allows to subtract non-flow effects from flow mea- 
surements order by order. Note that some experimental 
artifacts, such as track splitting, in the analysis also con- 
tribute to the two particle correlation; in this respect 
multi-particle techniques are also valuable, as they sup- 
press such contributions as well. 

One of the problems in using multi-particle correlations 
is the computing power needed to go over all possible par- 
ticle multiplets, which practically prohibits calculations 
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of correlations of order larger than k — 3 (three-particle 
correlations). To avoid this problem, it was suggested 
in [7] to express cumulants in terms of moments of the 
magnitude of the corresponding flow vector Q n , defined 
as: 
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where M is the number of particles. Unfortunately, flow 
estimates from cumulants constructed in such a way were 
systematically biased by the interference between various 
harmonics. An improved cumulant method using the for- 
malism of generating functions suggested in [8, 9 j fixed 
the problem of interfering harmonics while keeping the 
number of operations still linear with multiplicity M. For 
this approach the analytical calculations become rather 
tedious and therefore the solutions are obtained using in- 
terpolation formulae. Unfortunately this introduces nu- 
merical uncertainties and requires tuning of interpolat- 
ing parameters for different values of the flow harmon- 
ics v n and multiplicity. More recently a Lee- Yang-Zero's 
sum method [TTHT1] has been developed to suppress non- 
flow contribution to all orders. Closely related to that 
are methods of Fourier and Bessel transforms of the Q- 
distributions [TS] , and the method of direct fitting of the 
Q-distribution. All these methods, while indeed being 
almost insensitive to non-flow, are biased by interference 
of different harmonics. 

In this paper we present a new method to calculate 
multi-particle cumulants in terms of moments of (in gen- 
eral, different harmonics) Q-vectors. In our approach the 
cumulants are not biased by interference between various 
harmonics, interpolating formulas used in the formalism 
of generating functions are not needed, and, moreover, 
all detector effects can be disentangled from the flow es- 
timates in a single pass over the data at the level of or 
better than any other method. The number of operations 
required in our approach is still oc M for each k. Since 
in our approach cumulants are calculated without any 
approximation and directly from the data we often call 
them direct cumulants (also referred to as Q-cumulants 
because they are expressed analytically in terms of dif- 
ferent harmonic Q- vectors). 

Flow fluctuations are an important part of an 
anisotropic flow study. It is believed that flow fluctu- 
ations are mostly determined by initial geometry fluctu- 
ations 16J of the system created in a collision. An im- 
portant consequence of this is that the anisotropic flow 
develops relative to the so-called participant plane(s) in- 
stead of the reaction plane determined by the direction 
of the impact parameter [T7]. We note that the method 
to calculate cumulants proposed in this paper is not in- 
fluenced by how exactly the anisotropic flow is being de- 
veloped. 

In our simulations we show results obtained up to the 
8-th order cumulant, although we think that in practice 
there is little advantage to go higher than order six, be- 



cause going to higher order does not remove the system- 
atic uncertainty related to contribution from clusters ex- 
hibiting flow (see the discussion of systematic uncertain- 
ties associated with cumulant analysis in [IB])- For ex- 
ample, in a 4-particle correlation analysis this bias corre- 
sponds to the situation when two particles are correlated 
because they are coming from the same cluster and, in 
addition, correlated with another two particles via flow. 

The paper is organized as follows. After the main def- 
initions are introduced in section [TTJ we describe how the 
so-called reference flow can be calculated. The reference 
flow is an average flow in some momentum window; it is 
needed for the calculation of the differential flow of parti- 
cles of interest. To optimize the procedure, the reference 
flow can be calculated using weights, e.g. weighted with 
transverse momentum of the particle. Thus the reference 
flow can be noticeably different from integrated flow of 
the same particles. Section IV describes how the differ- 
ential flow is calculated. To show how the method works 
in different environments and how it compares to some 
other methods we show simulation results in section \V\ 
Finally, we summarize the main features of the method. 
Technical details, including the derivation of the main 
equations, equations in case of using non-unity weights 
in the calculation of reference flow, and acceptance effects 
are provided in Appendices. 



II. MULTI-PARTICLE AZIMUTHAL 
CORRELATIONS AND CUMULANTS 

In this paper we discuss mostly 2- and 4-particle 
azimuthal correlations (formulae for 6-particle correla- 
tion are provided in the Appendix), but the generaliza- 
tion to azimuthal correlations involving more particles is 
straightforward. The method can be easily applied for 
calculations of mixed harmonics multi-particle correla- 
tions. In fact, mixed harmonics correlations are needed 
in our approach for calculations of any multi-particle cor- 
relations with order higher than 2. Presenting 4-particle 
correlations below, we also show how the 3-particle cor- 
relations, involving one particle of a double harmonic can 
be calculated. All the correlations are obtained by first 
averaging over all particles in a given event and then av- 
eraging over all events. The latter may involve weights 
depending on event multiplicity. 

We define single-event average 2- and 4-particle az- 
imuthal correlations in the following way: 
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means that all indices in the sum must be taken different. 



3 



The second step involves averaging over all events: 



III. REFERENCE FLOW 
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where by double brackets we denote an average, first over 
all particles and then over all events. W/ 2 ) and Wu\ 
are the event weights, which are used to minimize the 
effect of multiplicity variations in the event sample on the 
estimates of 2- and 4-particle correlations. In general, 
the optimal choice of weights would be determined by 
the multiplicity dependence of v n . The best approach 
might be to calculate the cumulants at fixed M and then 
average over the entire event sample. In our calculations, 
with v n independent of multiplicity, we use: 

W {2) = M(M - I) , (9) 
W {4) =M(M-l)(M-2)(M-3). (10) 

The above choice for the event weights takes into account 
the number of different 2- and 4-particle combinations in 
an event with multiplicity M. 

The general formalism of cumulants was introduced 
into flow analysis by Ollitrault et al [7H!]- We will use 
below the notations from those papers. The 2 nd order 
cumulant, c n {2}, is simply an average of 2-particle cor- 
relation defined in Eq. Q: 



Cn{2} = ((2)) 



(11) 



As was pointed out first in [5] the genuine 4-particle cor- 
relation (i.e. 4-particle cumulant), is given by: 



c n {4} = «4»-2.«2)r 



(12) 



Expressions (II) and il2n are applicable only for detec- 



tors with uniform acceptance and will be generalized in 
Appendix [C] to extend their applicability for detectors 
with non-uniform acceptance. 

Different order cumulants provide independent esti- 
mates for the same reference harmonic v n . In particu- 
lar 0: 



v n {2} = ^c n {2} , 
v n {A}= ^-c„{4} ; 
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where the notation v n {2} is used to denote the reference 
flow v n estimated from the 2 nd order cumulant c„{2}, and 
v„{4} stands for the reference flow v n estimated from the 
4 th order cumulant c„{4}. 



To obtain the 2 nd order cumulant it suffices to separate 
diagonal and off-diagonal terms in |Q„| 2 : 

\Qn 



M 

1 e™^-^ = M+ ^ e m{ ^-^ ] , (15) 



which can be trivially solved to obtain (2): 
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The event averaging is being performed via Eq. ([T]) . The 
resulting expression for ((2)) is than used to estimate 2 nd 
order cumulant (see Eq. (11)), which in turn is used to 
estimate the reference flow harmonic v n by making use 



of Eq. @. 

To obtain the 4 th order cumulant we start with the 
decomposition of \Q n \ 4 (for details, see Appendix |a| 
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We have four distinct cases for the indices i, j, k and 
I: 1) they are all different (4-particle correlation), 2) 
three are different, 3) two are different or 4) they are 
all the same. Note, that the case of three different in- 
dices corresponds to the so-called mixed harmonics 3- 
particle correlations, in many analyses of great interest 
by itself |T5l [TS] . Equations for 3-particle correlations are 
provided in Appendix[Aj Taking everything into account, 
we obtain the following analytic result for the single-event 
average 4-particle correlation defined in Eq. 
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The reason why the originally proposed cumulant anal- 
ysis [7] was biased lies in the fact that the terms con- 
sisting of Q- vectors evaluated in different harmonics (for 
instance terms \Q2n\ 2 and [Q2nQ n Qn\) have been ne- 



glected. As seen from Eq. (18), such terms do appear in 



the analytic results and are crucial in disentangling the 
interference between harmonics. In particular, if a higher 
harmonic v 2n is present than |Q„| 4 picks up an addi- 
tional contribution depending on that harmonic, namely 
vl n M(M-l)+v n v 2 n2M(M-l){M-2), which is exactly 
canceled out with the contribution of harmonic v% n to 
\Q 2n \ 2 and 9te [Q2 n Q„Q n ], which read Muf„(M-l) and 
M(M-l)(M-2)vlv 2n +M(M-l)vi n , respectively. 

The final, event averaged 4-particle azimuthal correla- 
tion, ((4)), is then obtained by making use of Eqs. 
and (10). Using ((4)) and ((2)) one can calculate the 



order cumulant from Eq. (12) 
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The reference flow is mainly used to calculate differ- 
ential flow. Therefore, one can optimize the calculation 
of reference flow to minimize the uncertainties in the fi- 
nal results. This is done by using different weights (e.g. 
particle transverse momentum) in the definition of flow 
vectors used in reference flow calculations. We provide 
all the equations necessary for calculations with weights 
in Appendix |B) 

The equations so far are applicable for an analysis with 
a detector with full uniform azimuthal coverage. In a 
non-ideal case one needs to take into account the ac- 
ceptance corrections O US] ■ Acceptance affects the cu- 
mulants in three ways: (i) contributions from additional 
terms, e.g. proportional to ((cosn</>)) or ((s'mncf))), that 
for a detector with full uniform azimuthal coverage are 
identical to zero, (ii) contributions from other flow har- 
monics, and (iii) the cumulant might be rescaled, which 
at the end can affect the final extracted flow values. We 
refer to Refs. [T^J HO] for a more complete discussion of 
acceptance effects. In practice the most important cor- 
rection is the first one, for which we provide the full set 
of equations for a 2- and 4- particle cumulant analysis. 

The generalized 2 nd order cumulant which can also be 
used for detectors with non-uniform acceptance is: 

c„{2} = ((2)) - <Ke( [ «cosn0i» + i ((smnfa)) ] 



x [ ((cosn<^ 2 )) - i ((smnfa)} 



= ((2)) - ((cos,^)) 2 - ((sinn^)) 2 , (19) 

where for the last line we have used the fact that for in- 
stance ((cos n<f>i)) and ((cos rufa)) are the same quantities 
apart from the trivial relabelin g. R emarkably, only two 



additional terms appear in Eq. (19), namely ((cosn</>i)) 
and ((sinn</>i)) , which counterbalance the bias to ((2)) 
coming from very general detector inefficiencies. Fur- 
ther details on treating the acceptance effects, including 
formulae for the 4 th order cumulant are provided in Ap- 
pendix [Cl 



IV. DIFFERENTIAL FLOW 

Once the reference flow has been estimated with the 
help of the formalism from previous section, we proceed 
to the calculation of differential flow. For that, all par- 
ticles selected for flow analysis are labeled as Reference 
Flow Particle, RFP, and/or Particle Of Interest, POL 
These labels are needed because flow analysis is being 
performed in two steps. In the first step one estimates 
the reference flow by using only the RFPs, while in the 
second step we estimate the differential flow of POIs with 
respect to the reference flow of the RFPs obtained in the 
first step. 



A. Reduced multi-particle azimuthal correlations 

For reduced single-event average 2- and 4-particle az- 
imuthal correlations we use the following notations and 
definitions: 
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where m p is the total number of particles labeled as POI 
(some of which might have been also labeled addition- 
ally as RFP) , m q is the total number of particles labeled 
both as RFP and POI, M is the total number of par- 
ticles labeled as RFP (some of which might have been 
also labeled additionally as POI) in the event, ipi is the 
azimuthal angle of the i-th particle labeled as POI and 
taken from the phase window of interest (taken even if 
it was also additionally labeled as RFP), 4>j is the az- 
imuthal angle of the j-th particle labeled as RFP (taken 

even if it was also additionally labeled as POI). , as 
before, denotes the sum with all indices taken different. 

Final, event averaged reduced 2- and 4-particle corre- 
lations are given by: 
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In our calculations we use event weights w^') an d W(A') 
defined as: 
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B. Differential cumulants 
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We derive equations for the differential equations with 
the help of p- and g-vectors; the former built out of all 
POIs (m p in total), and the second only from POI labeled 
also as RFP (m q in total): 
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The g-vector is introduced here in order to subtract ef- 
fects of autocorrelations. Using those, we have obtained 
the following equations for the average reduced single- 
and all-event 2-particle correlations: 
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For detectors with uniform azimuthal acceptance the 
differential 2 nd order cumulant is given by 



dn{2} = ((2')) , 



(30) 



where, again we use notation from Ref. [5). We present 
equations for the case of detectors with non-uniform ac- 
ceptance in Appendix [C| 

Estimates of differential flow v' n are being denoted as 
w[j{2} and are given by [5]: 
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Below we present the corresponding formulae for re- 
duced 4-particle correlations: 
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The 4 order differential cumulant is given by [5]: 
d n {4} = «4'}}-2.«2'»«2» . 



(32) 
(33) 

(34) 



Equations for the case of detectors with non-uniform ac- 
ceptance are again presented in Appendix [C| 

Having obtained estimates for d n {4} and c n {4}, we can 
estimate differential flow [8]: 
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Similarly to reference flow, we use the notation i^{4} for 
differential flow harmonics v' n obtained from 4 th order 
cumulants. t | ' l {4} and v' n {2} are independent estimates 
for the same differential flow harmonic vL . 



V. SIMULATION RESULTS 

We have tested the new method with extensive sim- 
ulations. The results, presented below, show that the 
method effectively suppresses non-flow contributions, il- 
lustrate the ability to remove the interference of the dif- 
ferent harmonics, show the applicability for detectors 
having significant acceptance "holes" , and give an exam- 
ple of a differential flow analysis. In the figures, v 2 {MC}, 
shown in the first bin, represents the Monte Carlo esti- 
mate for v n , which was obtained using the known reaction 
plane event-by-event. Other estimates in the figures are 
obtained without using this information. 
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FIG. 2. Elliptic flow extracted by different methods for 10 5 
simulated events with multiplicity M = 500, V2 = 0.05 and 
at the same time 1)4 = 0.1. 



Figure [2] shows the results from a simulation of events 
with anisotropic flow present in two harmonics, the sec- 
ond and the fourth. Elliptic flow estimated by differ- 
ent methods is shown in the figure. A clear bias is ob- 
served in the estimates from fitting of the Q-distribution 
method and the Lee- Yang Zero's Sum method, labeled 
as u 2 {FQD} and w 2 {LYZS}, respectively. Results ob- 
tained with direct cumulants of different order, labeled 
as V2{k,QC}, are unaffected by interference. 

To demonstrate that the method works well even in 
cases with rather bad acceptance we simulated 10 7 events 
with v 2 = 0.05 for a detector that had two large "holes" 
(see Fig. [3^,) . Figure [3 o shows the obtained v 2 estimates 
using Eqs. (Ill and (121) which are valid for detectors 
with perfect acceptance using open markers. Clearly 
these values are strongly biased. The v 2 estimates ob- 
tained from the more general equations for cumulants, 
namely Eqs. (CI) and (C6), which do account for the 



acceptance effects are shown as closed markers and agree 
with the Monte Carlo estimate. In Fig. |3fc we look in 
more detail at the agreement with the Monte Carlo esti- 
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(c) 

FIG. 3. a) The azimuthal distribution of accepted particles, 
b) Extracted elliptic flow accounting for acceptance effects, 
closed markers, and without, open markers, c) Extracted el- 
liptic flow accounting for acceptance effects in different meth- 
ods. 



making use of particles labeled as RFPs (using Eqs. (Ill, 



(12), (13) and ( pL4]f) . The estimates of reference flow are 
presented in Fig. 141 In the second step we estimate the 
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FIG. 4. Reference flow extracted from particles labeled as 
RFPs (pions in Therminator) 



differential flow of POIs (in this example protons were 
labeled as POIs) with respect to the reference flow of 
RFPs estimated in the first step. For each p t bin we 

and 
flow 



evaluate d n {2} and cZ„{4}, and use equations (31 



( 35 1 to estimate differential flow. The differentia 



results for protons are presented in Fig. [5] The resulting 
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mate and, in addition, compare to other methods. 

The figure clearly shows that detector effects arc cor- 
rected for at the level of or better than other methods. 

As an example of a differential flow analysis we show 
results for v' 2 [pt) obtained with Therminator |21j . As 

RFPs we select pions and as POIs we select protons. p t -integrated flow of protons calculated by making use 



FIG. 5. Differential flow extracted for particles labeled as 
POIs from Therminator events (in this example we used pro- 
tons). The open circles denote 2 nd order estimate ( Eg. (31)) 
and closed squares denote 4 th order estimate (Eq. (35)). 
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of Eq. ([2]) is presented in Fig. [6j The figures for the 
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FIG. 6. pt-integrated flow calculated from Eq. Q of protons 
whose differential flow is presented in Fig. [5] 



integrated flow of the RFPs and POIs clearly show that 
the 2 nd order cumulant is biased by nonfiow while the 
higher order cumulants are in perfect agreement with the 
Monte Carlo. 



VI. SUMMARY 

In summary, we propose a new method to calculate 
multi-particle azimuthal correlations, which provides fast 
(in a single scan over the data) and exact (no approxi- 
mations) non-biased (no interference between different 
harmonics) estimates for cumulants. In the paper, we 
provide the corresponding formulae for correlations up 
to the 6-th order, but the method, if needed, can be gen- 
eralized for higher orders. 

We have not discussed issues of the cumulant approach 
in general, such as multiplicity fluctuations, flow fluctua- 
tions, and low sensitivity for small flow values, but believe 
that our method will be helpful in investigating all these 
questions. 

The proposed method has been extensively tested in 
simulations and has been used for real data analysis by 
the STAR and ALICE Collaborations [22l |23]. Fur- 
ther details about the method, including equations for 
8-particle correlations, equations for estimates and eval- 
uation of statistical errors, comparison to other methods, 
can be found in [2"3"l. 
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Appendix A: Equations for 3-, 4- and 6- particle 
correlations 



Below we use the following definitions: 

M 



tV,J , 



,i=l 
M 



M,2 



(Al) 
(A2) 



M 

in{24> i ~it> i -<t> k ) (A3) 



i,j,k=l 



W 2n\n,n ~ r> 7 , 

W/n.nl2n = W 
(^)n,n|n,n r> 



n,n\2n — \ %J /2n\n,n > 
M 



_ e in(<t>i+<j>j-<l>k-<t>i) 



^ ,J ^i.j,k,l — l 

Using this notation one finds: 

\Qnt = (4)„,„|„,„ ' PMA 

+ vV2n\n,n + \^/n,n\2n 

+ (2)„|„-4P M ,2(M-1) 

+ (2)2n|2n ' Pm.2 

+ 2P M ,2 + M . 



(A4) 
(A5) 



ai.3 



(A6) 



The 2-particle correlations (2) ra |„ was already expressed 
in terms of the Q- vector evaluated in harmonic n, see 
Eq. (fl6|): 



(2)2ra|2n — 



\Q 2 n\ ~M 



P. 



A/.2 



(A7) 



To obtain (3) 2jl | n n and (3) n n | 2n we have to decompose 

Q2nQ* n Q*„ = (3) 2n | n ,„--PM,3 + (2) n | n -2P W) 2 

+ (2) Pao + 1-M, (A8) 

and Q n QnQ*2n- After inserting resu lts for (2)„|„ and 
(2)2n|2rt given in Eqs. (16) and (A7), we arrive at the 
following equality: 



, = 2 ^ [Q2nQ* n Q* n ]-2-\Q n \ 

\ A In,n\2n T \' 3 /2n|n,n Z M{M - 1)(M - 2) 

|Q 2 „| 2 -2M 



M(M- 1)(M- 2) 



(A9) 
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After inserting Eqs. (16), (A7| and (A9) into Eq. (A6) 



n,n\n,n 



the 



and solving the resulting expression for (4) 
single-event average 4-particle correlations (Eq.(ll8j)) fol- 
lows. 

This derivation can be generalized to obtain analytic 
results for any higher order multi-particle azimuthal cor- 
relations. Below we provide the expression for the 6- 
particle correlation: 



<6> = 



1 



P 



MM 



M 



in{4>i,+<j>j+<f> h -<f>i-<l>m-<i>n) 



\Qn\ 6 +9 ■ \Q 2n \ 2 \Qn\ 2 ~6 ■ [Q 2 nQuQ* n Q* n Qn} 

M(M - 1)(M - 2)(M - 3)(M - 4)(M - 5) 

i 9te [QznQnQnQn] ~ 3 • Me [Q 3n Q$„Q* ] 

M(M - 1)(M - 2)(M - 3) (M - 4)(M - 5) 

9(M - 4) ■ 9U [Q 2 „Q T *Q;j + 2 ■ |Q 3 n| 2 
M(M - 1)(M - 2)(M - 3)(M - 4)(M - 5) 

O IQ»| 4 +IQ 2 n| 2 

M(M - 1)(M - 2)(M - 3)(M - 5) 
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\Qn\ 

M(M - 1)(M - 3)(M - 4) 
6 

" (M-l)(M-2)(M-3) ' 
With that, the 6 th order cumulant is given by 
c„{6} = ((6))-9-((2))((4)) + 12.((2)) 3 , 
and the reference flow v n is estimated as 



v„{6} 



r c„{6} 



(A10) 



(All) 



(A12) 



Appendix B: Particle weights 

Below we provide formulae to use for the case when 
the reference flow is calculated using particle weights. 
For that we introduce a weighted Q-vector evaluated in 
harmonic n: 



In.k 



M 



(Bl) 



where lUj is a particle weight of the i-th particle labeled 
as RFP and M is the total number of RFPs in an event. 
In general, we will need flow vectors with power k up 
to the order of multi-particle correlations. Similarly, we 
define 



Vn.k = ^2 



(B2) 



Note that only particles which have a RFP label, have 
a non-unit weight, while for the particles labeled as POI 



only, Wi = 1. For the subset of POIs which consists of 
all particles labeled both as POI and RFP (m q in total) 
we introduce 



q n ,k = E 



For RFPs we also introduce: 



' M 

E* 

.i=l 

M 



Mabcd - = wfw^wf ■ 

i.j 7 k,l,... — l 



(B3) 

(B4) 
(B5) 



For all particles labeled both as RFP and POI we evaluate 
the following quantity: 



E 



W: 



(B6) 



while in the definition below the first sum runs over all 
POIs in the window of interest and the remaining sums 
run over all RPs in an event 



m D M 



M' 



bed-- 



J2 y! w >H w ?--- 

i=l j,k,l,..—l 



(B7) 



Using the definitions presented above the weighted 
single-event 2- and 4-particle correlations are given by: 



M 



WiWj e 



in(<f>i — <f>j ) 



M 



w.WjWkWi e 



(B8) 

in(4>i+4>j~4> k -4>,) (B9) 



ij,k,l—l 

The event weights ^ and ( 10 ) now read 



W {2) = Mu , 
= Mim 



(BIO) 
(BH) 



Analogously, the reduced single-event multi-particle cor- 
relations now read: 



m v M 



( 2 ')^EE'-v 

■^oi i=1 j=1 



in(ipi — <l>j) 



A I 



1 I 

^^m 1 — E E w ^ w/ ' 

0111 i=lj,k,l = l 



(B12) 

in{tl>i-\<t>j-4>k-<l>i) (B13) 



where the event weights ( 24 1 and ( 25 ) are now 



w {2 , } = M' m , 
w< 4 ') = Km 



(B14) 
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The weighted average 2-particle correlations are given by 
the following equations: 



(2) 
((2)) 

Mu 



\Q 



n,l\ 



Si, 2 



^2,1 — Si, 2 

EliiMnh 

M 

- ' r, "'j 

: 5*2,1 — Si, 2 , 



(B15) 



and the weighted average 4-particle correlations are given 
by: 



(4) = 



IQn.il 4 + \Q2nM 2 - 2 • sac [q 2 „, 2 q;,iQ;,i] 
s-mc [Qn,zQ* n ,i] - 4-Si, 2 |g„,i| 2 



— 6 ■ Si, 4 — 2 • 5*2,2 
M 

•Mini = ^2 w l w J w k wi 



I Mi 



((4)) 



i,j,k,l=l 

■ Si t i — 6 • S± t 232,l + 8 • Si,3Si,i + 3 • S2,2 

6 • S 1A , 

E"i(Mim),(4), 
E^i(A4iin)i ' 



(B16) 



where the weighted Q-vector, Q n ,fc, was defined in Eq. 
plj) and S p , k in Eq. p3|. 



Weighted reduced 2- and 4-particle azimuthal correla- 
tions are given by the following formulas: 



(2') 
«2')> 



PnfiQLl 



}pSi, 



Sl,l 



Er=i(^'oi)i 

X' 01 = ^ W 3 : = m pSi,i - Si,i 

2—1 2, J — 1 



and, 



(4') = PnfiQn,lQn,lQn,l 

— Q.2n,lQ*n,lQ*n,l ~ PnfiQn, lQln, 2 

- 2 • Si t 2Pn,oQn,l ~ 2 ' Si,i |Q„,l| 2 
+ 7 • q n ,lQn,l ~ Qn,\1 n ,2 

+ l2n,lQ*2n.2 + 2 • PnfiQ n ,3 



2 ■ Sl,lSl 2 — 6 • Si 



/■Mom , 



EtiKm). (40 
ELCM'oiii). 

771*, M 



«4'))= 



-Mom=5Z W J W k w l 
i=l j,k,l=l 

=m p [S 3 ,i - 3 • Si,iSi, 2 



2 • Si, 3 ] 



-3-[si I i(S 2 ,i-S 1) 2)+2-(s 1 , 3 -si,2Si, 1 )] (B18) 

We note that to evaluate all quantities appearing on the 
right hand sides in analytic expressions (B15 B18 1 only 
a single loop over data is required. 



Appendix C: Non-uniform acceptance 

Building cumulants from multi-particle correlations we 
have so far omitted terms which vanish for the detectors 
with uniform acceptance. For a more general case they 
have to be kept [7J [5J US \M ■ The more general 2 nd order 
cumulant now reads: 



Cn{2} = ((2)) 



((cosn</>i)) 2 + ((sinnc^i))' 



(Cf) 



The correction terms can be expressed in terms of the 
real and imaginary parts of the Q-vector Q: 



Eli (^He [Q n ])i 



((cosrc^i)) = 



((sin ncpi)) 



Eti Pm[gn])j 

EL^ 



(C2) 
(C3) 



When particle weights are used the average 2-particle 
correlation ((2)) is determined from Eqs. (B15), while 
Eqs. ([C2l and (C3) generalize into: 



((cosn^i)) = 



((sinn</>i)) = 



Eti (^[Q»,i]) t 

Ei=i(Si,i)i 

Eti Wft4 
Ei=i(Si,i)» 



(C4) 
(C5) 



(B17) where Q„ i can be determined from the definition of the 



weighted Q-vector (Bl| and Si,i from definition (B4| 
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The generalized 4 order cumulant reads: 
Cn {4}H(4))-2-((2)) 2 - 

- 4 • ((cOSTK/li)) ((cOS?l(0i - 02 - 03))) 

+ 4 • ((sinn^i)) «sinn(0i - 2 - <fe))) 

- ((cosn(0i + <fe))) 2 - ((sinn(0i + 2 ))) 2 
+ 4- ((cosn(<?!>i + 2 ))) 

x ((cosn</>i)) 2 — ((sinn^i)) 2 

+ 8 • ((sinn(^i + <j> 2 ))) ((sinn^i)) ((cosra^i)) 
+ 8-((cosn(<?!>i-0 2 ))) 



((cosn^i)) +((sinn0i)) 
6- ((cosn^i)) 2 + ((sinn^i)) 2 



(C6) 



The terms starting from the second line in Eq. ( C6 ) 
counter balance the bias coming from non-uniform ac- 
ceptance so that c n {4} is unbiased. These terms can be 
expressed in terms of Q- vec tors: 



// r± r j. \\\ Ei=l (St* [QnQn Q2n\)i tr<n\ 

«coB»(*x+fc)»= ^ {Mi _ x) • (C7) 



((sinn(0i+0 2 )))= 



Ei=l ( Jm [Q«Qn-Q2n]). t 

E^iAf^Mi-l) 



(C8) 



and the Eqs. ( C9 1 and (CIO) generalize into 



JV 



((cosn(0 1 -0 2 -0 3 )»=j Y, [Q n ,lQ*n,lQn,l 
^ i=l 

1 - 

-Qn,iQ2„,2-2-5i, 2 Q;, 1 +2-Q; 3 ]). /^(MuOi 

J t=l 

f N 

((sinn((/>i-^2-^3)))=-| ^ (3m [Qn.lQn.lQn.l 

-Qn,iQk 2 -2^i, 2 Ql 1 +2-Q;, 3 ]) i /^(M m )i 



M 



Mih=5ju; < «; j -ii;a^S , 3 ) i-3-Si,25i,i+2 -Si )3 . (C12) 

i,j,k=X 



The generalized 2 nd order differential cumulant reads 



JV 



^ i=l 
1 N 

-2(M-l)5te[Q^) < /^Ii(I r l)(Mr2), (C9) 

J z=l 

r w 

^ »=1 

-2(M-lptn[Q;]) i /^M,(M r l)(M,-2). (CIO) 



d„{2} = ((2'» - 
((cosra/'i)) ((cosn</> 2 )) — ((sinm/>i)) ((sinn</> 2 )) (C13) 



Expressions for ((cosn^i)) and ((sinn^i)) were already 
given in Eqs. (C2| and ( |C3[ ), respectively (when particle 
weights are being used in Eqs. (C4| and (C5|, respec- 
tively). Similarly: 



When particle weights are used the average 2-particle 
correlation ((2)) is determined from Eqs. (B15), the av- 
erage 4-particle c orre lation ((4 )) is determined from Eqs. 
(|B16|), the Eqs. |C7| and ([C8t generalize into: 



Ei=l (^ e [Qn,lQn,l — <32n,2] 



((cosn(0i + (/> 2 ))) 
((sinn(0i + (/> 2 ))) 



Ei=l (^ m [QnsQnS ~ Q 2 n, 2 ]) j; 



EL(^u) 



,1/ 



((cos nipi)) 
((sin nipi)) 



Eili(mp). 
EL (Jmb ra ]) t 

Eili(w p )i 



(C14) 
(C15) 



where p n and m p were denned in Section IV The Eqs 



.M11 = tUjWj = S2 1 — 5*1 2 , (Cll) (C14) and (C15I remain unchanged when particle weights 



: j=i 



are being used. 
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The generalized 4 order differential cumulant reads: 
d n {4} = ((4'))-2.((2'))((2)) (C16) 

- ((COS 77/01)) ((COS 71(01-02-03))) 

+ (<sinn0i)) ((sin7i(0i-0 2 -0 3 ))) 

- ((cOS770i)) ((COS 71(^1-02-03))) 

+ ((sinri0i)) ((sin n^i- 02 -0a))) 

- 2 • ((cOS7J0l)) ((cosn(?/>i+02-03))) 

- 2 • ((sinn0i)) ((sinn(0>i+0 2 -03))) 

- ((cOS77(0i+0 2 ))) ((cOS7l(0i+0 2 ))) 

- ((shi77(0i+0 2 ))) ((sinn(0i+0 2 ))) 

+ 2- ((COS 77(01 +02))) 

x [((cosn0i)) ((cos7i0i)) — ((sinm/>i)) ((sinn0i))] 
+ 2- ((sinn(0i+0 2 ))) 

x [((cosn0i)) ( (sin 7701 )) + ( (sin n^i)) ((cosn0i))] 

+ 4- ((COS 77(01 -02))) 

x [((cosn0i)) ((cos ti0i)) + ((sinn0i)) ((sinn0i))] 

+ 2- ((008 77(701+02))) 

x ( (cos ?70i )) 2 — ((sin 7701 )) 
+ 4 • ((sinn('0i+02))) ((cosn0i)) ((sin 7700) 

+ 4 • ((COS 77(01 -02))) ((cOS?70l)) 2 + ((shl770l)) 2 

- 6- ((cos 770i)) 2 — ((sinn0i)) 2 

x [((costt/0i)) ((cost70i)) — ((sinn0i)) ((sinn0i))] 

- 12 • ((cos?70i)) ((sinn0i)) 

x [((sinn^i)) ((cos 770i)) + ((cos 770i)) ((sin?70i))] 



r 

((008 77(0^02 -03))) = E [Pn^Q^-PnQan] 

^ i=l 

-9te [2m 3 Q;-2^] ). [(m p M - 2m g )(M - 1)]. , 

' i=l 

((sin 77(01 -02 -03))) = \ 2_j {^ m [PnQ*nQn - PnQ* 2n } 



N 



-3m{2m q Q* n -2q* n }). \/Y,[{m P M~2m q ){M-l)] i . 
i=i 

(C19) 



When particle weights are used Eqs. (C17) generalize 
into: 



((COS 77(01 +02))) 



((sin 77(01 + 02))) = 



11 , k - Q2n,k])i 



(C20) 



The terms starting from the second line in Eq. (C16) 
counter balance the bias coming from non-uniform accep- 
tance. Some of the new terms appearing in this expres- . . 

sion can be expressed again in products of flow vectors: < CIS ) generalize into. 



I^lgte \p n Qn-q2n]) i 



((COS 77(0i+0 2 ))) = 



((sin 77(0i+0 2 )))= 



(C17) 



((COS 77(01+02 -03))) = \Y,( mC [Pn(\Qn\ 2 -M^j 

^ i=l 

-SKt [q2nQn+m q Q n -2q n ] ) . W ^ [(m p M -2m q )(M - 1)]. 



i=l 



((sin 77(01 + 02 -03 



))) = E(3m[p„(|Q„| 2 -M 
^ i=l 



((COS 77(01+02 -0 3 ))H X! (lQ«^| 2 - 5 

-$He [g2n,lQn,l + s l,l<3n,l - 2 <7n,2] 

JV 

^ (m p (5 2 ,i - ^1,2) - 2 ■ (31,151,1 - si )2 )) 

i=l 

((sin7l(0i+0 2 -0 3 ))H (Jtn (|Qn,l| 2 - S 'i i 2) 



-3m [^n.lQn.l + Sl.l^n.l -2 ?"^] 
-3m[92nQ;+m g Q n -2 gn ]).^/^[(mpM-2m,)(M-l)] i , w 

(018) I j=1 J 
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and finally, Eqs. (C19) generalize into 



JV 



((cOSn^l-02 -0 8 ))H E [Pn (Qn,lQ;,l-^2„, 2 )] 

-2- We kiQn,i - «;,a] )<}/ 

j E ( m p( S 2,l - S l,2) - 2 ' ( S 1,1 5 1.1 ~ S l,2))j f i 
^ i=l J 

r N 

^ i=l 

-2 -3m [si.iQ^i-^]),}/ 

[ £ K(5 2 ,i - Si l2 ) - 2 • (si,iS M - sid). ] . (C22) 
^ i=l J 
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